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_ ' An analytic solution of the Hartree-Fock problem for a 2DEG at filling 1/3 and half an electron 

per unit cell is presented. The Coulomb interaction dynamically breaks the first Landau level in 
• three narrow sub-bands, one of which is fully occupied and the other empty, as in the composite 

fermion model. The localized orbitals associated to the Bloch like single electron wavefunctions 
are nearly static, resembling the angular momentum eigenstates within a Landau level for non- 
interacting fermions. Strong correlations are expected owing to the large charge density overlap 
between neighboring plaquettes. A numerical evaluation brings the cohesive energy close to that 
of the best present day models. It is also found that correlations are long range, requiring over 
50 particles spread over a finite sample to approach convergence. Since presently allowed exact 
calculations are far from this number, the question of how relevant the considered wave-function is 
for the description of the ground state of the 2DEG system remains open. 
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It is well accepted today that the fractional quantum Hall effect results from the formation of a strongly correlated 
spin polarized electron liquid in the lowest Landau level, that crystallizes below filling l/ 7b 2 i 3 . It is intriguing, 
however, that a mean field charge density wave solution exists at all fillings, called a Hall crystal (HC), which yields 
in a natural way the odd-denominator rule that characterizes the effect 4£ The state is characterized by a unit cell 
that contains a fractional charge of even denominator, such as e/2. The detailed investigation of this class of mean 
field states was stopped largely due to their higher energy than that of the liquid, and the very small correlation 
correction obtained by Yoshioka and Lee (YL) for the Wigner solid (WS), which accomodates one whole electron per 
plaquette. 6 However, while in the latter case there is little overlap of electrons from neighboring cells, in the former 
the charge density forms continuous ridges between cells opening the way for an appreciable increase of the cohesive 
correlation energy^ This idea was first explored in the simplest case of 1/2 filling £ The results confirmed the effect 
yielding a second order correlation correction an order of magnitude larger than for the best WS state. The same line 
of thinking was also recently considered by Mikhailov.(— ). 

In this work we report results for the correlation energy of the HC at 1/3 filling. The work rests on an analytic 
solution of the mean field problem constructde using symmetry considerations. We show that the electron-electron 
interaction breaks the single particle degeneracy of the lowest Landau level (LLL), splitting it into three remarkably 
flat bands. All bands contain the same number of states so that only one is filled while the other two are empty, 
separated by a large gap. While the HC unit cell traps 3/2 flux quanta of the original field, it is only pierced by one 
flux quanta of the reduced field B* = B 1 / 3 — B 1 / 2 , where B 1 / 3 and B 1 / 2 are the fields at 1/3 and 1/2 filling factor, 
respectively. The narrowness of the occupied band and this latter fact suggest that the magnetic Wannier states are 
approximate solutions of the problem and behave as nearly free quasiparticlcs filling the LLL of the effective field B*, 
much as the composite fermion model predictsAiS 

We compute the correlation energy using the YL method, that is, evaluating the second order perturbation correction 
to the Hartree-Fock (HF) solution. Two methods are used, one employing a Monte Carlo approach for the space 
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integrals over a large sample, and the other by computing in a discrete momentum space using a set defined by 
periodic boundary conditions applied to a comparatively small sample. As for the half filling case, both methods 
yield an energy correction of about 10%, an order of magnitude larger than that for the WS stated, confirming the 
earlier conjecture that at all fillings the high electron overlap of the HC enhances significantly the cohesive energy 
corrections 7 '-, thus re-opening the case for this state as the proper precursor in a first principles perturbative approach 
to the true gound state of the system in the thermodynamic limit. 

One further outcome of our calculations is that over 50 particles would be required in order to approach convergence, 
way above current permissible sizes used in first principles numerical evaluations normally performed with up to 12 
particles, number that has been extended to treat at most 27 particles^. This may be understood again looking at 
the large inter cell overlap through hexagonal ridges, which seems to point out the importance of cooperative ring 
exchange effects in the system^ To be operative, these need large sample sizes and would be missed by such numerical 
diagonalizations of the hamiltonian. If the true ground state retains the essential features of its HF approximation, 
then one would expect the Hall crystal to be the proper mean field precursor due to its higher correlation energy than 
other mean field solutions^ 

In Sec. 2 the single particle Hartree-Fock hamiltonian is diagonalized analytically. Section 3 is devoted to the 
second order correction os the mean field solution. Finally, in Section 4 we present our conclusions. Details of the 
calculation are given in appendices A-D. In particular, in Appendix B a formula is derived expressing the single 
particle Hartree-Fock Hamiltonian in terms of the magnetic translations associated to an arbitrary periodicity lattice 
of the HF problem. This formula, up to our knowledge, is absent in the literature. 



II. BLOCH REDUCTION OF THE HF PROBLEM AT v = | 

A. One particle Hamiltonian and magnetic translations 

We consider a 2DEG in a strong perpendicular magnetic field. We are interested in an analytic treatment of the 
Hartree-Fock problem at filling 1/3. We start out by writing the Hartree-Fock hamiltonian in the compact form 

ff H *. = 5>(Q)exp(-^)T r 2 nxQ , (1) 

Q 

where the coefficients w(Q) are given by 

„(Q) = 2^p(Q)exp(^) r^ex P (-^) - ^I (^)) A (2) 



4 ' \r Q |Q| ±v 4 ' V 2 v 4 7 e Q r, 

where e Q is the background dielectric constant. The magnetic translation operators T are defined in Appendix A, 
while form (1) is derived in Appendix B. The particle density in real space p(x) is assumed to be periodic under 
displacements covering the triangular lattice defined by the vectors 

R = r 1 a 1 +r 2 a 2 , n,r 2 = 0,±1,±2, (3) 

/ 6-7T . 1 \/3\ 

a2 = V71 ( 2'- )ro - 
The Fourier components of the density are defined as 

P(Q) = [ dx p(x) exp(zQ.x), (4) 

Acell J 

where A ce u is the unit cell area 

A ceU = n.aixa 2 = 3^. (5) 
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Here n is a unit vector normal to the electron gas plane. Through this cell traverse a flux | in units of the magnetic 
flux quantum <j) — hc/e. Associated with the above real space lattice is the reciprocal lattice 
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B. Block diagonalization of the HF hamiltonian 

We next show that it is possible to find a basis in which the matrix representation of (1) has a diagonal form 
composed of simple 3-dimensional blocks. In addition, the functions have such a structure that they automatically 
furnish the translation symmetry of the total density over the lattice (2). For this purpose we consider the basis 
functions </? P (x) defined in Appendix A, constructed over the lattice with primitive vectors 

bi = ai, (7) 

, 2 

b 2 = ^a 2 . 

The magnetic translations then have the form 

T r l nxQ = r -|Q 2 bi+Qib 2 - ( 8 ) 

Since the flux piercing the unit cell of the lattice (2) is not an integral number of flux quanta, the set of translation 
operators obtained by varying the integers Qiand Q 2 in © do not commute and it is not possible to find common 
eigenfunctions to all of them. The basis we shall construct decomposes in a set of three dimensional subspaces, closing 
each of them under the action of translations for all values of Q. 

A first step in finding the basis is to define a set of eigenfunctions of a translation in the vector a 2 , which is a period 
of the density. For a given value of the momentum p such orbitals may be written as 

Xp (x) = -L (V P (x) + rA rTa^p(x)) , a = ±1. (9) 

p V2 V cxp(-za 2 .p) J 

Using the formulas in Appendix A it can be readily proven that these functions satisfy the eigenvalue relation 

T^x^x) = aexp(-ia 2 .p) *£(x). (10) 

The two values of a appearing in these equations will play an important role in what follows. They will allow us to 
impose the periodicity of the density under the shifts in ai and a 2 , in spite of the impossibility of obtaining a basis of 
eigenvectors of all the magnetic translations in the lattice. It should be stressed that the range of values of p defining 
independent functions in the new basis have been reduced in half, the two values of a compensating for this reduction. 
The restriction comes from the singular property of the basis {</? p } , that a magnetic translation is fully equivalent to 
a shifting of the momenta argument as implied by the relation 

Tr^p(x) = •7 7 p(R.)¥'p+ttA(R)( x )> (H) 

where JF p (R) is a pure phase factor. See Appendix A and Ref (16) for the justification of this relation. It directly 
implies that the shift done in a 2 in constructing the new basis precisely changes the momentum p of the particle in 
— Si/2. Therefore, what has been basically done is to form linear combinations of the original functions associated 
with different values of the momentum. In the new basis the magnetic translation in ai has the simple effect 

T aiXp "(x) = exp(-iai.p) Xp CT W, (12) 

that is, it merely changes the sign of a and multiplies it by a phase factor. 

As the next step let us employ the fact that, although the functions are not eigenvectors of translations in a l5 
they are eigenfunctions of the double sized translations in 2ai. This is because its effect, when considered as two 
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consecutive shifts in ai, have the simple result of making two consecutive changes of sign of a that reproduce the 
original function. Therefore, if for p and a fixed we construct the triplet of states formed by the function Xp( x ) 
and the other two obtained by a pair of successive translations in the vector -2ai/3, the operation of performing a 
translation in an arbitrary multiple of these vectors leaves the triplets invariant. 

A specific linear combination within each triplet which is also an eigenfunction of the translation in -2ai/3 can be 
obtained by constructing the new basis 



Ap 



(x)= Yl c r(p) T -| sbl Xp(x), r = -1,0,1. 



(13) 



After solving the linear equations for the constants c s r obtained by imposing the condition that these functions be 
solutions of the eigenvalue problem 



T_ 



)1 X P r ' CT) W = Axp riff) (x) 



one finds 



A< r >(p) = exp(^p.a 1 + ^) 

c = 7f exp ^3 ip ' aiS ~ 

r,s = -1,0,1. 



(14) 



2irirs 



Substituting in (fT*5f) yields 



X p r ' CT) (x) 



r 
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P 



-7= J] exp(--ip.a 1 s — )T_2 sbi X p(x), 



-1,0,1 



±1 
P + 



s i , s 2 n j-1 xo 

n— + m 7p n i m = 0, ±1, ±2, . 



(15) 



(16) 



The last relation expresses that, modulo a phase factor, the states of the new basis are equivalent upon a shift of 
their momenta p in any linear combination with integer coefficients of half the unit cell vectors of the reciprocal 
lattice corresponding to the periodicity of the density. The periodicity of the states under the shifts in si/2 was 
discussed above, and the one related with S2/2 similarly follows from the relation (|53|l in Appendix A, expressing the 
equivalence of a magnetic translation with a shift in momentum. The functions just defined have an alternative and 
more compact form given by 



* 5>xp(iP(P 



r,cr) 



5ni 



m 1 m 2 )T m (j)(-x.) 



A, 



(3,2) 



e = £ 1 (sc 1 ) + e 2 (2c 2 ), 



(17) 



(18) 
(19) 



where the effective momenta and the new elementary lattice of vectors m over which the sum is performed are given 
by 



a- 13 
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2 2' 



(20) 



miCi + m 2 C2i 



1 3 ' 
a 2 

C2 = — • 
3 

The double sum Ijl7|l can be evaluated to obtain an explicit formula in terms of the Elliptic Theta functions as follows 
cx P (-^4) 



Xp r ^(x) 



(3,2) 
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where 
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))9)|ri), (21) 



Let us inspect now the action of a magnetic translation by |a 2 on the new functions. If such a transformation 
leaves the triplets invariant, then the matrix reduction of the Hartree-Fock Hamiltonian will follow. One has, 



T |a2 x P r ^(x) = -^ Y, exp(-^p.a lS -^)T ia2 r_ |sbiX -(x). 
But after using (|49() for changing the order of the two operators within the sum, it follows that 



(22) 



r ia2 X P r ' c7) (x) = A p (-a 2 )x p 



([r-l],*) 



(x), 



(23) 



where the square bracket defines the number among the set {-1,0,1} that is equivalent, modulo 3, to the integer in 
the argument. 

For fixed p and a the matrix elements of the Hamiltonian Q in the new basis can readily be found to have the 
form 



= J2 «(Q) exp(-^) exp | / po , Q r + !—(),[,-■ 

Q 



2m 

T 



Qi 



V,[r-Qi]- 



(24) 



The problem has thus been reduced to the self-consistent diagonalization of a three dimensional matrix for each value 
of momenta p and parameter a. The basis can be checked to have the following set of transformations properties, 



(r,cr 



T 2c2 X {r ' a 

J 3ci Ap 

^3C2 Xp ' 

ix 



(r,(7 



)(x 

)(x 
)(x 
)(x 
)(x 



cxp(-2z p . Cl -^p) Xp ^)(x), 

exp(-2* p.c 2 ) Xp [r - 11,CT) (x), 
ex P (-3* p.cjv^- )(x), 
erexp(-3i p.c 2 )x p r ' <T) (x), 

(— r,cr) / \ 

X-p ; (x), 



(25) 



where the parity transformation / is defined as usual by / Xp'^i^) — Xp' a \~ x ). As shown in Appendix A, from 
these symmetry properties it follows that the density associated with any Slater determinant constructed by selecting 
one orbital within each triplet has exact periodicity under shifts in all vectors R. 

In order to find the solution of the mean field problem by an iterative technique it is sufficient to make an ansatz 
for the density in the first step, and then diagonalize numerically the matrices for a sufficiently high partition of 
the reduced Brillouin cell momenta p defined by (| 1 fc>|) - By selecting the normalized lowest energy state within each 
three dimensional quantum mechanical problem, the Fourier components of the density corresponding to the new step 
should be constructed. Following its definition, it can be done by means of the following expression, 



P(Q) = ^~ E /dx|5> (p) Xp ^ ) (x)| 2 exp( l Q.x) 

p,cr— ±1 r 

VS]>i ~~"E E e*p(ip.nxQr2 + ^Q a Q 1 )x 



(26) 
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Qa)£r ',[r+Qi]- 
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This formula can be obtained by evaluating the Gaussian integrals appearing after substituting the expansions defining 
the functions Xp • The coefficients g®, g\ and g%, define the components of the eigenvectors of the single particle 
HF Hamiltonian in the basis of states \- They fix the wave-functions of the filled band and the empty bands as 

^ M) (x) = <?r(p)x£ ,CT) (x) (27) 

r=-l,0,l 

6 = 0,1,2. 

Here 6 = labels the filled band in each triplet and 6=1,2 label the two empty bands. 



C. Functions that vanish at the origin are exact solutions 



In order to proceed within an analytical context let us consider the observation from former numerical studies, that 
the particle density when the number of flux quanta piercing the unit cell is a half integer essentially vanishes at all 
lattice points^ Then, let us first assume that the density rigorously vanishes at this set of points. If so is the case, 
the wave-function of any of the filled states should then also vanish at those points. This requirement follows from 
the fact that the Hartree-Fock particle density is a sum over the individual densities of all occupied orbitals, 



£ lE-^pk^wi 2 ^*)- (28) 

p.cr— ±1 r 

We can then use this property to fix the coefficients of the wave functions within each triplet. 

After imposing the vanishing conditions at the origin, the coefficients defining the functions l|27|) are fully determined 
and take the form 



30°(P) 



A/p* 



x P °' +) (o) ^'(oj-xp^fo) xJr w) (o) 



y - l(V> Kf* (-l,+l)/ n N (+-l)/r.\ (+l,+l)/ n N (-1,-1) / n s [ J) 

;v pXp (0) X P (O)-Xp (0) X P '(0) 



i xp <- 1 \o)xp- h+1 \o)-xp°> +1) (o)xp- 1 <- 1 \o) 
^ x { P h+1 \o) X p +h - 1} (0) - x ( P +h+1) (0) xp-'^io) 

1 = I.9o°(p)| 2 + I.9°i(p)| 2 + I3i(p)| 2 . 

Note that the coefficients g are all independent of a. This completes the definition of our functions. That they are 
true solutions of the Hartree-Fock problem is proven in Appendix C. 

The particle density may now be computed replacing these functions in Eq. 1)28(1. The real space particle density 
thus obtained is shown in Fig^ 

A main property to be noticed in this figure is the formation of sharp hexagonal channels surrounding the low density 
regions at the center of which the vanishing density occurs. These structures mark the difference with the Wigner 
solid whose unit cell encloses one flux quanta. The charge density in this latter instance is made up essentially of well 
localized gaussians centered at each lattice point. In our case there is strong overlap, suggesting that cooperative ring 
exchange involving many unit cells is a large contribution to the correlation energy. 

Further, the insertion of the calculated density in the eigenvalue equation associated with the matrix representation 
of the Hamiltonian in each triplet, Eqs. an( i C2J, allows for the evaluation of the one particle spectrum of the 
system. As it was expected, three energy bands appear, each associated with a value of the index r and covering the 
full range of p within each triplet. We also note that states associated with a = ±1 turn out to be degenerate. The 
bands dispersion relations are illustrated in Fig. [3 

Note the narrowness of the bands as compared with the gaps separating them. This fact leads to the idea that in 
this mean field approximation the Coulomb interaction, although breaking the degeneracy of the first Landau level, it 
reorganizes the states in three equally populated separate sets that again are aproximately degenerate, as if they were 
Landau levels of a renormalized problem. Since at 1/3 filling one band is full and the other two empty, one expects 
the electrons to behave dynamically similarly to filling one, except for a different effective mass, as the composite 
fermion model suggests»2^£ The magnetic Wannier states are expected to be approximate solutions, playing the role 
of the angular momentum states in the non-interacting problem at filling one. 



FIG. 1: The density of particles associated with the HC state. Its structure indicates the high overlapping of the localized 
Wannier-like states in terms of which the HF determinant always can be expressed. 
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FIG. 2: Triplet of narrow bands in which the Landau level is splitted by the action of the Coulomb interaction. 

These properties seem to corroborate the possibility of tracing a link between weakly interacting composite fermions 
and the Bloch or Wannier orbitals in the mean field solutions considered earlier by one of us4£ To finish this section 
we note that the mean field energy per particle obtained from the above formalism for the solution we have constructed 
confirms the value reported in an earlier numerical computation^ 



-0.362- 



(30) 
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III. CORRELATION ENERGY IN SECOND ORDER 

In order to obtain the energy correction to second order of perturbation theory we proceed similarly as in Ref (8), 
starting with the expression^ 

£ P) = J2^ HF I ( h - h hf) I ^) e hf_ e . ^ I ( h - h hf) I * ffF ) (31) 

i 

v I (* HF j_g_j_gj) I 2 

i 

Here , £'- ff - F are the mean field Slater determinant and total Hartree-Fock energy associated with the ground 
state, respectively, and H is the projection of the exact many particle Hamiltonian onto the first Landau level. The 
many particle excited states VE^ are Slater determinants constructed with the basis (|27|l . mixing states in the filled 
band with states in the empty bands. It follows that | = 0, a property that allowed to write the last 

equality in (|31|l . In the second quantized representation the Hamiltonian H will have non- vanishing matrix elements 
linking the HF state and excited states of the form | $,•) = a v ey a^, \ Q HF ) only where creates an electron 
of quantum numbers £, etc. The index i is a shorthand notation for the set of two pairs of filled (77, rf S F) and 
empty (£, ('€T-f) electron states, where T and F are the set of all states in the Landau level and the filled ones, 
respectively. The indices 77 =(0,p,cr) and £ =(b,p,a) for b = 1,2 denote the quantum numbers of filled and empty 
states, respectively. The total energy of an excited state is Ei = E^ HF ^> + e(£) + e(£') — e(rj) — e(r/). Then, the second 
order correction can be rewritten in the form& 



F(2) = v v i(^|g w 44i^> I 2 m] 



where the total projected Hamiltonian is 

2e a 



H = / / <W ^( X )**(x')^-^*(x')*(x) (33) 



° a,a' /3, 0' 

The matrix elements of the Coulomb interaction are given by 

M(a,a' \P',f3) = J J rfxdx' $* (x)$*, (x') r^^ri 'M'O^M, (34) 

where we have used the shorthand notation $ Q = $p' <T ' ) . By using the anti-commutation relations [o a ,o^,] = 6 a , a >, 
formula l|31|l can be expressed as 

(2) = v v \ii dxdx 'n.v'(^)j^^'^')\ 2 



where the two particle states <& v ,r)' are defined by 

<I>, /( x)<T> mxOj-_<IK(x><Ivx') 
The pairs (77,7/)and (£,£') are considered as unordered. 



*W (x, x') = ?n " ^ " ■ (36) 



A. Correlation energy: first evaluation 



As pointed out above the single particle bands are remarkably flat. We can use this property to simplify the 
calculation of the energy correction. First, we approximate the filled band energies appearing in the denominator of 
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(|35|l by their mean value, 

In addition, and slightly more crudely, we substitute the energies in the excited bands by a common energy equal to 
half the sum of the mean energies of the two bands, 

e(0 = —j— 

This last approximation is taken in view of the small relative gap separating these two bands. 
With this simplifying substitution (|35|) can be expressed in the simpler form 



g//" A // W |x-x-||y-y-| X (37) 



2(2e _e 1 _e 2 )e2 
(t/ (x' , y' ) TT f (x, y ) - TT f (x' , y ) 717 (x, y' ) ) 7r e (y , x)7r e (y' , x' ) , 

where iTf is the projection operator on the subspace of states of the filled band and 7r e the projector associated to the 
subspace of states formed with the union of the empty bands. The projectors have the following expression in terms 
of the corresponding densities 

tt/^x') = ]>>(^)(x)($(°^(x')r, (38) 

p,<T 

= P(x, x'^vrr 2 . £ p f (Q) exp(»Q.(^±iL + ^nx(x- x'))), 
Q 

7T e (x,x') = ^0>^)(X)($(^)(X'))*, 

P,fa,(7 



F(x, x'^vrr 2 £ p t (Q) ex p(iQ.(^-^- + \i n x (x - x'))), 



Q 

1 

27rr; 



^ e (Q) = 7^7- 2^Q,o — Pb(Q), 



where ^q,o is the ordinary Kronecker delta and P(x, x') is the projection operator onto the first Landau level, defined 
in Appendix A. These expressions can be obtained from formula <|71|> in Appendix B. After evaluating a few spatial 
integrals, the following formula for the correlation energy per particle e' 2 -* = /N is obtained 



e (2) = 9 , 9 ^ , 2 47r ^ r ° (7rr ° 2)2 EEEP/(QiW(Q 2 )^(Q3)Pe(-Q 1 -Q 2 -Q3) (39) 
Z(Zeo-ei-~e2)e v — " — — 



Qi Q2 Q 3 

eM iQl+ f )2r h J dz-l=exp(-^)exp(-(Q 1 +Q3).^-z(Q 2 +Q^. 2 
I - n x (Qi+CbK - z(Q 2 +Q 3 )r ) 2 

(exp i(— - n x (Qj+Qgjro - i(Q 2 +Q 3 H) 2 - exp(™(— - n x (Q x +Q 3 )r - i(Q 2 +Q 3 )r ) 2 )). 

O To o T 

Further progress in the evaluation of the energy correction is aided by noting the presence of an exponential factor in 
the squared Fourier wave- vectors, allowing to take just a few of them only, for sufficient convergence. We thus assume 
the 36 Fourier components associated with the shortest values of Q are nonzero, only. Of these, those associated with 
the longest Q are slightly altered to assure that the sum rule^ 

Q 

is satisfied. This condition is used to cancel various fictitious divergences in the formula defining the correlation 
energy^. Then, we obtain 
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FIG. 3: Second order correction to the energy per particle as a function of inverse number of sample points in the Monte-Carlo 
integration. 



£(2) = T-#7^^EEE A /(Qi) A /(Q2)A e (Q 3 )A e (Q 1 + Q 2 + Q 3 )x (41) 
v 1 2 > Qi Q 2 Q 3 



r »,„ , ^ , _ fr% , ^ x _ , (Q 2 +Q 3 ) 2 ^, 



exp(-i 7^(Q 2 + Q3) • n x (Q ± + Q 3 )r 



2 



x 



, / (u-i(Q 2 + Q 3 )r Q ) 2 . (u-i(Q 2 + Q 3 )r ) 
du exp z r (Q 2 +Q 3 ). 



7 ( g (u - i(Q 2 +Q 3 )r ) 2 ) sinh ( -(u - i(Q 2 +Q 3 )r ) 



^Q2+Q 3 ,o 



'(u + n x (Qi + Q 3 )r ) 2 



where the order parameters of the filled and empty bands have been defined as usual, 

A,(Q) = 27rr^ / (Q)cxp(i^), 



A e (Q) = 2^ Pe (Q)exp( 
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The integrals in Eq. 141(1 were calculated using the Monte-Carlo method with a variable number of sample points 
in the range N — 10000 — 40000. The region of integration was taken as a square of size L = 50 in units of the 
magnetic length ro, centered at the origin. At 1/3 filling this sample contains 132 particles and traps 396 flux quanta. 
The results for the correlation energy as a function of 1/N are shown in Fig. El In units e 2 /e r the estimated value 
for the correlation energy is 

e (2) ~ -0.04 ±0.01 . (42) 

Adding this correction to the Hartree-Fock energy l|30|) one obtains the value e ~ — 0.40± 0.01, to be compared with 
the result of Yoshioka and Lee with its own correction e ~ —0.394, and that of Laughlin e ~ —0.416, all in the same 
units. 



B. Correlation energy: dependence on the sample size 



The lack of a definite convergence pattern provided by the method just described prompted us to try a different 
procedure. Using the formulas derived in Appendix D it is possible to recast Eq. 1411 to give it the form, 
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^ = jM- E E E E ( 4 3) 

1 

y 

e (0, Pfj, op) _|_ £ (0, pp, a p i) _ e (b, p a , cr Q ) _ e (b', p a ,, a a ,) 
| ^ K \V^i, r0l +Pp, rfi -Pa,r a -Pa',r a ,,0) X 

^(P a ,fa.^) 5*k(Pa)^*(Pc«',?V,0v) ^(Pa/^P^'.r^,^/) 5^, (P/3' ) X 

^(pp^trp) g° rp (jp )^„ Jn x (P^-Pa.rjr^) x 
^^(nx (P r ^,-P„', v )r2) x V(a, a', /?,/?') | 2 , 

where the function V is defined by 

V(a, a', /?,/?') = E^T15 4? 7^exp(-r 2 (P^,-P QVa ,+Q*) 2 )x 



exp [-i r 2 n x Q* ■ (P a , r<s + P^ - P^, r(J) -P a /, r<i( ) 



The form of the pure phase factors J 7 * and J-p , the momenta P Q . r and the special reciprocal lattice vectors Q* are 
all specified in Appendix D. As before, the functions g%(p) are the coefficients determining the single particle HF 
excitations. We remark that once the coefficients for the filled band are known the other set of coefficients for the 
empty bands can be evaluated analytically as the two eigenvectors of the (3x3) matrix representation of the single 
particle HF hamiltonian, orthogonal to the vector <?°(p). 

In order to evaluate numerically expression l|43|) we restrict the Hilbert space of the single particle HF problem 
by defining a cell of sides 2JViai, 2N 2 & 2 with N\, N% integers, and imposing periodic boundary conditions over its 
borders. Using the symmetry properties (|25|l 



these constraints can be translated into the relations 



Tn^X^'X*) = exp(-z6p.c 1 7V 1 )x P r ' ff) (x), 
^ 2 6c 2 X^ CT) (x) = exp(-z6p.c 2 7V 2 ) X ^ ) (x). 



-6p. Cl = 2^—, 

n 2 

-6 P .C 2 = 27T— . 

They restrict the values of the quasimomentum p to a discrete set, as expressed in compact form by the condition 

n L t 1 _n 1 t 1 

v Ni 2 N 2 3 v ' 

< Tlx < N\ ,0<n 2 < N 2 , 

where ti and t 2 are defined in Appendix A. The cell chosen traps 6NiN 2 flux quanta and contains 2N\N 2 particles. 
The integrals over the continuum in the Brillouin zone of Eq. (|43[1 are now restricted to summations over a finite 
number of values of p. Formula 1)43(1 was first evaluated by for an elongated sample with N\ = 1, and N 2 = N 
varying from 1 to 13. In order to keep the computing time within our allowed practical limits, the values of n\ and 
n 2 in the expression defining the summation argument 

3 2 

Q* = m- ti +n 2 - t 2 , (45) 
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FIG. 4: Correlation energy as a function of the parameter N. The circles indicate the calculated values for elongated samples 
of the form (Ni, N2) = (1, JV) . The black dots signal the values corresponding more appropriate samples with form defined by 
(Ni,N 2 ) = (N,N). 



were restricted to the regions —4 < n\ < 4 and —4 < n.2 < 4. The results of this calculation are shown as empty circles 
in Fig. 0. Note that the correlation energy per particle grows with N, approaching rather slowly the thermodynamic 
limit whose value e^ 2 -* ~ —0.06 e 2 /e r exceeds the result ()42(l of the previous method. The energy correction was also 
calculated for samples with an aspect ratio close to one, taking (N 1: N 2 ) — (N, N), up to N = 5, an upper bound set 
by computing limitations. The results are shown in Fig. by filled circles. As can be observed, the energy values 
for those samples grow similarly with N than those associated with elongated samples. This property indicates that 
these more appropriate samples also need values of N greater than five in order to approach convergence. However, 
already at N — 5 the number of particles involved is 2N 2 = 50, considerably higher than the numbers handled in 
exact diagonalization. 

The most relevant finding of this last method of evaluation is the need for very high sample sizes in order to 
approach the thermodynamical limit. This outcome strongly supports the relevance of long distance correlations in 
the attainment of the true character of the ground state through numerical computations. 



IV. SUMMARY 



An analytic solution of the Hartree-Fock problem for a 2DEG at filling 1/3 and 1/2 electron per unit cell is found by 
employing symmetry considerations and a special complete set of common eigenfunctions of the magnetic translations. 
A triplet of bands arises that turn out surprisingly flat as functions of the two dimensional quasimomentum. The 
Coulomb interaction breaks the first Landau level in three narrow sub-bands resembling effective Landau levels of 
composite fermions, one of which is fully occupied and the other two, empty. The energy per particle reproduces 
an earlier numerical result for this quantity. We find that the charge density forms hexagonal rings throughout the 
lattice, suggesting strong long range correlations owing to cooperative ring exchange. An evaluation of the correlation 
energy using second order perturbation theory yields a correction an order of magnitude larger than that for the 
localized single-particle features of the usual Wigner Crystal stated and lowers the energy to make it comparable with 
the lowest values obtained using other models. Our results also suggest that the thermodynamic limit in a cell with 
aspect ratio near one is not reached if less than 50 particles are included. 
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VI. APPENDIX A 
A. Eigenfunctions of magnetic translations 

Consider a 2DEG constrained to move in a plane of area A under a perpendicular magnetic field B. A useful basis 
set of single particle Bloch-like states in the lowest Landau level can be defined in terms of linear combinations of the 
normalized zero angular momentum eigenfunction 

«*-lk-«-V!> (46) 

in the compact formii^*iSi ? , 




N k =^N 4>0 jy J-l)^exp(ik.£. _ 



.2 ■ 

a 



Due to its role in the above definition the function <f> is called the "seed" function. The sum runs over all integers 
£±,£2 defining a planar lattice L, through £ — £±hi + £2^2, where the unit cell intercepts one flux quantum, so that 
n.bixb2 = 2irrg. The magnetic translation operators Tr acting on any function / are defined by 

2?V 

T R /(x)=exp(— A(R).x) /(x-R), (48) 
he 

where the vector potential is assumed in the axial gauge A(x) — B(—X2,Xi,0)/2 and the electron charge e is taken 
with its negative sign. This basis was employed before to obtain exact mean field solutions of the related problem at 
filling 1/2& For arbitrary vectors Ri and R2 the translation operators satisfy the commutation relation 

T Rl T R2 = exp(-A(R 1 ).R 2 )T R2 T Rl . (49) 
As it may be easily verified, the functions (pk satisfy the eigenvalue equation 

T t ^ k (x) = A k (£) ^ k (x), (50) 
AkO?) = {-lf l12 exp(-z k.£). (51) 

Arranged in an arbitrary Slater determinant these functions are exact solutions of the Hartree-Fock problemiiSiii 
This strong property happens because the HF single particle Hamiltonian commutes with all translations leaving L 
invariant^ The functions (|47fl are common eigenfunctions of the commuting magnetic translations. Moreover, the set 
of eigenvalues l|51|l uniquely determines them. Therefore, the HF hamiltonian associated with the Slater determinant 
can not change those eigenvalues and the y?k should be eigenfunctions. 

Finally, let us show that the effect of an arbitrary translation on the basis functions is equivalent to a shift in 
the momentum label, modulo a phase factor^ Operating twice with the translation operator involving an arbitrary 
vector a and a vector in the lattice £ and using Eqs. I)49|l and H5(J[) one readily gets, 



T a T e (^p(x) = Ap(^)I> p (x) (52) 

Again using H50JI we have 



exp(2 — A(a).£) T e T a ^ p (x). 



TiT a v? p (x) = A p+27feA(a) (^)T a <p p (x). 
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Then, taking into account that the set of eigenvalues defines uniquely the wave-functions modulo a phase, it follows 
that 

T a </J p (x) = Jp(a) <yS p+27 fLA(a)(x), (53) 

from which also follows, 

^pW = 77T- ( 54 ) 

^p+2 7ft :A(a)(a) 

That is, a magnetic translation is equivalent to a shift in the quasi- momentum. 

B. A ^-transform 

Any function / in the first Landau level, and its inverse, can be represented as 

pes 



dp 

c(p) tfpW, 



c(p) = J rfx ^*(x) /(x) 



where <£> p (x) = y/S~ip p (x), c(p) = Vd c(p), with 5 = I-ktq iV0 o , with A^ the number of flux quanta in the system 
area, B is the Brillouin zone defined by the unit cell vectors 

tl = J 11 X b 2, 

t 2 = -^n x bi. 



The orthogonality and completeness relations in the first Landau level of the modified functions take the forms 

/ dx$£(x)$y(x) = (2^) 2 <5(p-p'), 

P(x,x') = ]T ^pW ^(x 7 ) = / 7^ ^p(x) ^ p (x') (55) 
pes J #* 

C. The orbitals xp as special cases of the basis functions <^ P (x) 

Lets us now verify that the functions Xp constructed to diagonalize in blocks the HF Hamiltonian in Sec. 2 are 
simply constant phase factors multiplied by a particular kind of basis functions ip p . These differ from the orbitals 
considered in Sec. 2 and defined by the unit cell vectors JJJ) in that their basis vectors are given by 

bi = -ar, (56) 
b 2 = a 2 . 

To start out it will be useful to consider the ip p -transform of a translated "seed" function 4>(x). For this purpose we 
notice that it is possible to fix the second argument of the projector operator to be equal to the translation vector a. 
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It follows that the projector operator in the first Landau level can be rewritten as 



P(x,a) = -l^exp(-^ i -^)exp(^A(a).x) (57) 

= J-^ r a0(x) = f £ p (x) <p*(a), 

\/2nrl J p eB \2^Y 

where the last equality is merely the definition of the projector operator in coordinate space. Therefore, these relations 
determine the following (^ p -transform of the translated "seed" function 

dp 



T a 0(x) = J ^ ^ ^27rr^;(a)J $5 p (x). (58) 
Next we will consider the definition 

and represent the vectors m in terms of alternative vectors £ and the indices a and f3 defined through the expressions 

m = 1 + 8, 
I = i x 2 ci +i 2 3 c 2 , 

5 = U Cl + U C2 . 

The symbols a, /?, ^ and ^ 2 are specified by the relations 

_ , x r if m 1 ={) Modulo{2) 
u - (mi) — 1 j j/ mi = i Modulo{2) ' 

-1 £/ m 2 = -1 Modulo{3) 

v = [mi] = { if m 2 = Modulo(3) , 

1 £/ m 2 = 1 Modulo(3) 

mi — ( m i) 

2 ' 
to 2 - [m 2 ] 



*i = 

^2 = 



3 

The use of these alternative definitions allows to write (|59() in the form 

X p r ' CT) 0r0 = ^ ( 3, 2 ) E(~ 1 )' 1 ' 2 expftP^.l) T t x (60) 
I E E exp(^ W + iP(^> CT ).(u Cl + vc 2 ))T UCl+VC2 <f>(x) 1 . 

Now, it can be noticed that the above expression differs form the definition (|47|l only in that the "seed" function is 
changed by a superposition of itself translated to the points itci + t>c 2 . However, it was shown in Ref. (16) that the 
basis functions are completely independent of an arbitrary change in the "seed" for any other function in the first 
Landau level. 

An alternative proof of this curious property is possible. To see it let us represent the " seed" functions in (|60|l by 
their </> p -transform (|58|l . yielding 



Xp* CT) (*) = r (32) E E ex P (—uv + iP^.(u Cl + vc 2 ))x 
/ (2$P V^^ 1 + ^)^« E exp(i(P<™"> - 
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after employing the character of eigenfunctions of the operators Ti that the functions ipq have. But now the following 
identity can be employed 

^exp(i((q-q')i)=^5(q,q') ) 



eeA 



where the 5 function is non- vanishing for equal arguments modulo a vector of the reciprocal lattice Q = §Qiti + |Q 2 t2- 
This relation allows to show the desired connection between the basis functions x^'^and </? p (x), 

where the phase factor is given by 

F{P,r,cr)= V V exp(^ww + iP (p ' r ' CT) .(uci +wc 2 )) (^p ( p, r , CT) (uci +i;c 2 ) . 



D. Periodicity of the Slater determinants of triplet orbitals 



Here we show that the single particle density associated with a Slater determinant formed with a set of functions, 
each one corresponding to an arbitrary linear combinations of the triplet of orbitals characterized by the indices (p, a), 
is periodic over the lattice R. To consider this question, let us use the basic property that the density of a many 
electron Slater determinant constructed with orthogonal orbitals is the sum of the individual densities of each orbital. 
Writing the selected linear combinations as 

*(r)(x)=^c(^) x (r^(x), (6i) 

r 

where the coefficients C p r are arbitary, and considering the symmetry properties it follows that 



„(x) = E E i^^wi 2 

p a=±l 

= E E K CT) (x)i 2 

p a=±l 

p a— ±1 r r' 

= E E EE C p r ' Cr) * C P^ ) ( ex P(£ A ( a i)- x )^p'^ ) ( x - a i)) exp(-|A(a 1 ).x) X (:^(x-a 1 ) 

p a— ±1 r r' 

= E E I^Hx-aOl 2 

P <T=±1 

= p(x-ai). 

One can show in a similar way that p(x) = p(x — a 2 ). 



VII. APPENDIX B 



Here we shall derive the general formula expressing the single particle Hartree-Fock hamiltonian as a sum of 
magnetic translation operations. The arguments of these operators are spatial vectors determined by the reciprocal 
lattice associated with the Fourier components of the density. It is worth remarking here that the lattice considered 
in this Appendix is not restricted in any way. Therefore, the representation of the single particle HF hamiltonian is 
valid for a general though periodic mean field problem. 
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Consider the action of the kernels of the direct and exchange interactions defining the HF hamiltonian over any of 
the elements (p p of the complete basis functions as follows 

^ (x) = H<i\ M y') = & y fdy dy> F(x, y')^-^^ P (y'), (62) 

j i y y i 

t ] (x) = H%\ ^ p (y') = -- fdy dy> P(x, y') p(y ' y p (y'), (63) 

e»i I y - y I 

where no is the jellium background charge density making the overall system neutral. The representation we shall 
derive is valid for both direct and exchange kernels and will be discussed below for each case separately. 

A. Direct Coulomb interaction 

The potential term representing the direct Coulomb interaction can be written in its Fourier transform representa- 
tion as follows 

„M(x) = ±f & ^ x ')-"°/ e2 



So J | x - X' | 

o prf n J I x - x I 



_ E ^£(QK^a2) exp(iQ . x) 

Q 
Q 

Then, using the following representation of the kernel for the projection operator in the first Landau level 

P(x, x ) = ^ exp( exp(-A(x ).x) 



we obtain, 

~ (d) 

(x) = j ay nx,yj 

Q 

Recalling definition (|47|l we obtain for each translation term 



(x)=y dy' P(x,y')E v(d) (Q)exp( J Q.y')^(y')- (65) 



4 d) (x) = f dy' P(x,y')J2v {d) (Q)eMiQy)Ti 4>(y'), 

J Q 

L_ £ W («0(Q) exp(-lr 2 Q 2 ) cxp (i^) exp ( - * (x - (* - r^Q x n)) 2 ), 
27rr q 1 2 4r 



lX .(£ + r^nx Q)), 

= E w(d) (Q)cxp(-ir 2 Q 2 )exp( l ^)T £+r g n)<Q ^(x). (66) 
Q 

This last relation can be further simplified by employing relation 152(1 and the similar property 

ic 

T Rl T R2 = exp(— A(R 1 ).R 2 )T Rl+R2 , 

he 



leading to the formula 
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f dy> P(x,y')^« (d) (Q)exp( J Q.y')r^(x) - ^ W W(Q) cxp(-ir 2 Q 2 )T r , nxQ T^(x). 
Q Q 

After adding corresponding the integrals over all translations t one obtains 

Hf F ^p(x) = ^ W ( d )(Q)cxp(-ir 2 Q 2 )r r g nxQ ^ p (x). 
Q 

Since the equality is valid for any element of the complete basis <p p (x) it then follows, 

H ( H d) F = ^ (d) (Q)exp(-lr 2 Q 2 )T rgnxQ . (67) 
Q 

B. Exchange interaction 

The derivation of the analogous representation for the exchange interaction kernel is more involved and needs for 
some special properties of the one particle density matrix. Then, let us initially consider these properties. 

1. One particle density matrix transformations 

The definition of the one-particle density matrix in terms of the Slater determinant $(xi, x 2 , x 3 , x 4 , ..~x N ) can be 
transformed by performing a simultaneous magnetic translation operation in a vector R of the periodic lattice over 
all the particle coordinates. Since this is assumed to be a symmetry transformation of the system, this map should 
leave the many particle state invariant. Thus by assumption 

P( x i> x 2) = J -J dx 3 dx 3 ..dx N |$(xi,x 2 ,x 3 ,x 4 , ..x w )| 2 (68) 
= / ... / dx 3 dx 3 ..dx N | J|T R (x i )$(xi,x 2 ,x 3 ,x 4 ,..x JV )| 2 

•* •* i=3 

= ]TvMxi)^(x 2 ), 

= Y,Tr(^M(T r (x 2 )^Jx 2 ))* , 

= ^exp(^A(R).( Xl - x 2 )) exp(*^(xi - R) (^„{x 2 - R))* p(xi,x 2 ), 

= exp(— A(R).(x 1 - x 2 ))p(xi - R, x 2 - R), 

an expression that furnishes the transformation law of the one particle density matrix under spatial shifts in the 
vectors of the periodic lattice. 

2. The density determines the density matrix 

Below we will show that the Fourier components of the density completely determine the whole one-particle density 

matrix. For this purpose, let us use the new variable z = x — x' in the one-particle density matrix, so that p(x, x') =p 
(x, z), where the tilde over p underlines the different functional expression of the new definition. It now transforms 
under spatial shifts as 

P (x - R, z) = exp(-— A(R).z) P (x, z). (69) 
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Therefore, the function 

g (x,z) = exp(-— A(x).z) P (x,z), 

is fully periodic in the variable x under the lattice shifts. That is 9 (x — R, z) =<? (x, z). Then, Fourier expanding g 
leads to the following expression for P 

P (x, z)=J2 cxpWQ - j- A(R))-x) p (Q, z), 
Q ° 

where 

p (Q,z) = -— / dx'exp(-jQ.x')exp(i^-A(z).x') P (x',z). (70) 
^4 ce /; J nc 

Further, let us consider the density written as the sum 

(*0=£ £ / nCT (p) X ^ CT) (x)^' j(T ,(p)x^ V) (x')*, 

and the fact that each x p r ' <T ' ) in the lower Landau level can be expanded as a linear combination of the complete basis 
functions {Ti </>(x)} defined over a lattice with one flux quantum per unit cellii^ Then, a generic term in the integral 
in l|70l ) has the form 

plj; (Q,z) = -— / dxexp(-iQ.x)exp(i-^-A(z).x) T t 0(x)(l>0(x') | X '= X _ Z )*. 

Since the function is a Gaussian, its magnetic translations H48|l have also this character, and thus the integrals in 
(170(1 can be calculated explicitly. We obtain, 

z 2 1 

pw (Q' z ) =pw (Q i °) cx p(-7t ~ TrQ-( n x z+i z ))- 

But since the phase factor here is independent of the indices £, £' , the summation of all the contributions in (|70() leads 
to the relation 

p(Q,z) = p(Q)exp(-^2 - -LQ.(n x z+i z)), (71) 
p (Q,0) = p(Q), 

which expresses the interesting result that the full one-particle density matrix is completely determined by the particle 
density. 

We now consider the exchange term l(r]5)l . Using the above results the integrals in 

(x) = -e 2 J dy' P(x, y') J dz £ cxp(z(Q - ^A(z)).y')^|^ ¥ , p (y' - z), (72) 

can be calculated explicitly. For this purpose, let us consider again the integral of a generic term Tg <p in the sum 
over £ defining the functions (p p . The integral for this term is then reduced to two simple integrals, which after some 
algebra can be explicitly evaluated to obtain 

i( £ )(x) = ^« (e) (Q)exp(-ir 2 Q 2 )exp( ? ^)r £+r 2 n><Q 0(x), 
Q 

= ^ U ( e )(Q)exp(-ir 2 Q 2 )r r g nxQ r^(x), (73) 
Q 



V (Q) = iTV 2 P (Q,0) exp( 4 r ° Q } /o( I r ° Q } - (74) 
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Henceforth, after adding all terms for different values of £ follows 



Hg F <p p ( X )= ^« (e) (Q)ex P (-ir2Q 2 )T r g nxQ j ^(x). 



Using the independence and completeness of the basis formed by the (p p 's in the first Landau level, we thus obtain 

H { h ] f = £^ e) (Q)exp(-ir 2 Q 2 )T r , nxp . 
Q 

Combining this results and the previous one (|67|l the following representation for the one-particle Hartree-Fock 
Hamiltonian in the first Landau level follows, 

Hhf = J2( vid) (Q)+v {e) (Q))oM-\rlQ 2 )T r , nxQ 
Q 

= ]>>(Q)exp(-ir 2 Q 2 )T r 2 nxQ , (75) 
Q 

with the coefficients defined by 

v(Q) = 27 rr 2 p(Q) ( ^ffi - exp(ir 2 Q 2 ) I (^ 2 Q 2 )) (76) 
VIII. APPENDIX C 

In this section we prove that the eigenfunctions that vanish at all latice points R are solutions of the HF equations. 
For this purpose, let us consider the Fourier components 

Py(Q) = -T— I dx p y (x)exp(-iQ.x) 

Acell J 

of the configuration space periodically extended density operator 

Py (x) = <$(x-y)=£«5(x-y + R) 

R 

= -r— y2 exp(-iQ.y) exp(iQ.x). 

Acell ~ 

We will examine the projection of this operator in the first Landau level. For this aim it is useful to consider the set 
of operators 

Pi = di + i-^-Aiix) 
he 

z ' 

g 

Gi = di- i—Ai{x) 
Tic 



^ 0^2 e V X 
z '0 



-1 

1 

where the indices i,j take the values 1,2. These quantities satisfy the commutation relations 
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[Pi,Gj] = 0, 

[Pi,Pj] = ~~J"> 
r O 

With their use, the spatial Fourier modes can be identically represented as 

1 



exp(-iQ.x) = cxp(-|(P - Gi).e ik Q k (2rl)) 
1 



= exp(--(P-G).nxQ(2r^)), (77) 

where the 3D-vectors P and G are defined as (Pi,P2,0) and (Gi,G2,0). But the two components Pi are linear 
combinations of the rising and lowering operators which transform functions of a Landau level in functions of the 
higher and lower contiguous levels. Therefore, the projection onto the first Landau level of the Fourier component can 
be obtained by simply omitting the operator P in the above formula. Then, the projection of the density operator 
will be 



p(f)(x) = I i-^exp(»Q.y)exp(-iG.nxQ(2r 2 )) 
Q 

= ~T~ ^2 eX P(*Q-y)^nxQ(r=) 
A ceU Q 



X* u 

where a new set of spatial lattice vectors was defined in terms of the Fourier wavevectors, 

x*= n x Qr-Q. 

Let us consider now the identity 

Po$>°(p) |X^ CT) >=0, 

r 

where the three coefficients g® (p), r = 1, 2, 3, are fixed by the condition that the linear combination of the functions 
of the two triplets of functions Ixp^' 7 ') for a = ±1 vanish at the origin x = 0. Projecting the above identity onto the 
first Landau level by acting with P at the left it follows 

Ppo ]T g° r (p) | 4^) = PpoPP J2 9°Ap) I 

r r 

= P { o P) T,9r(P)\4' a) }=0, (78) 

r 

implying that the projected periodic density has the defined functions as eigenvectors with zero eigenvalues. Let us 
inspect further the commutation properties of the projected density operator. Using the representation of the single 
particle Hamiltonian in terms of the magnetic translations found in Appendix B and the relations 

T x ,T y . = exp(£A(x*).y*)T x , +y ,, 

1 P 

T y ,T x , = cxp(--A(x*).y*)T x , +y ,, 



it follows that 
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Hhf,P^\ = £$>(x*)[T x .,T y ,] (79) 

x* y* 

= ^^ U (x*)(exp(^A(x*).z*)-cxp(-£A(x*).z*))T z ,. 

Considering that u is the coefficient of the superposition of translations giving the single particle HF Hamiltonian 
which is a function of Q 2 , and using the inversion symmetry of the solution under consideration u(x*) = it(-x*), one 
finds 



Hhf,P 



0. (80) 



Hence, acting with the Hamiltonian on (|78() and employing the definition of the projection operators and the periodic 
density operator it follows that 

HhfpW ]T <? r (p) | xt a) ) = PoHhf 9riP) I 

r r 

= f d a; 'P(x,x')po(x')^^(p)x^ ) (x') 

= f &'P(x,x')^«5(x'-R)^/i r (p)x^ CT Hx') 

J R r 

= ^P(x,R)^^(p) x ^)(R) 

R r 

= E ^3 cx P(-^f^) exp(£A(R).x) ^ Mpk^W 

R. r 

R V 27rr Q r 

= 0. 

Finally taking into account that the set of functions Tr<P(x) for all values of the lattice J2J are linearly independent, 
it follows that 

EMp)Xp' CT) ( R ) = °- (81) 

r 

In addition, since the coefficients of the linear combinations of the triplet functions for the two signs of a are uniquely 
defined by the vanishing condition of the function at the origin (or any point of the lattice R), we have, 

h r (p) = eg°(p), 

where e is here a normalization constant. Finally we obtain 

H HF J29r(p) I 4' a) )=J2 h r(P)\X^), (82) 



= zJ2gr(p) I Xp 



That is, the special functions that vanish at all laticc points arc eigenf unctions of the one particle HF Hamiltonian. 



IX. APPENDIX D 



In this Appendix we sketch the derivation of the formula i|43|) employed in the evaluation of the sample size 
dependence of the correction to the energy per particle. As before, we use the shorthand notation 
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r 

d (K,r,* a ) = T{pairj<7a) g K {pa) 

P a>r = p(P^) =pa _ rSl + ^_li3 S2 

The relations expressing the functions Xpa^^ m terms of <p Po times a phase factor has also been employed in order to 

define the new coefficients dp"' r ' a °'\ By representing the Coulomb potential by its Fourier transform, and substituting 
all the expansions of the functions ip Pa in terms of magnetic translations of the "seed" in the matrix elements (13411 . the 
resulting Gaussian integrals for each term can be evaluated. After some algebra it is possible to obtain the expresion 



M(a,a' | P',P)= I /dxd^$;(x)$* / ( X ') T ^-^ rT * /J /(x / )$ /S (x), (83) 



l 



dxdx' <^, Qr Jx) (Pp a ,^ (x') ; - _ — ; y PgVy (x') ip Pprj3 (x), 

E E E E (<: ,ra ) 45' :r °' ^ } ) * d S' r ,3 ' ' ff3 ' } ^ ,CTfl } x 



<S(P0',r„/ +P/3,r /3 -Pa,r«-P«', V ,0) X 

^> X (Py, V -P°',r a >o) ^ Q >(~ n X ( P y, V - P a ',v)^) 

/vr(2,3) ^(2,3) 

£ |P s ,,,,-pL. + Q-| + Q- » 2 ) 

exp [-i rgn X Q* • (P a ,r a + P/3,r/3 - P/S'.^-Po'.r,,,) 

where the summation vectors Q* and the norm functions are defined as follows, 

3 

Q* = m- si +n 2 s 2 , (84) 



<' 3) = v / ^jE(- 1 )^ 3ex p( iPJ -^) 



I = £ 1 {2c 1 )+e 2 (3c 2 ), 



and, as before, 5(P' , P) equals one when the arguments arc same modulo vectors of the class Q* defined above, and 
zero otherwise. It is possible to show now that the ratio multiplying this latter function a simply a phase factor. For 
this purpose let us consider the definition of the y>p , 

from which follows directly a connection between the value of this function at the origin and the norm functions 
appearing in (|83|l 



PP..r„ (°) = t/E ^ CXP ^ P^-^-^) = 



FT ^ 



Further, let us recall the formula connecting a translated function of the basis ip p and the same function for a shifted 
momenta (I53|l . After considering this relation evaluated at x = 0, the following expression can be obtained 
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^„.( nX (%,- p «'.v)'o) = ¥>P Q ,. Q ( n x (P a , ro -P fl , r/ ,)r§) 

= -^Pc.jnx (P^, r? -Pa,r a )?-o)Vp airct -||A(nx(P c «, ra -P )3(ry3 )rg)( )- 

But the new momenta argument of the function in the rhs can be simplified as 



to produce 



"l^( nX ( P a,r„-P/V>o) = ^ f X (n X (P^ -P^ )r Q 2 ) 



^P„„ Q (n x (Pf,, >r -P a >,r a ,)t$)=rp a , r Jn x (P^-Pa.rjrg)^ (0). 



Similarly it also follows 

^^ v (" nx ^.'^^'v)^) 3 ^,,^^ (V- p «'.v)'o)^, v (°) 
which finally leads to the desired results 

^P-.r a ( n x ( P /3,r,- p «^)^o)^P a ,. rQ ,(n x (P^^-Pa'^Jrg) 

The expression for the matrix element (|83|l allows to directly write formula l|43|> employed in Section 3 for evaluating 
the correlation energy. 
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